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ABSTRACT 


A  goal  of  The  Observing  System  Research  and  Predictability  Experiment  (THORPEX) 
Pacific  Asian  Regional  Campaign  (T-PARC)  that  was  conducted  in  2008  over  the 
western  North  Pacific  was  to  examine  the  targeted  observation  for  improving  initial 
conditions  of  numerical  forecasts  of  typhoon  track  and  intensity.  Results  indicated  that 
observation  impacts  depended  on  the  location  of  the  observation  with  respect  to  the  storm 
and  the  altitude  from  which  the  observation  provided  a  profile  of  winds,  temperature,  and 
moisture. 

In  this  thesis,  the  experiment  results  were  augmented  by  testing  the  impact  of 
observations  that  could  potentially  be  obtained  from  an  unmanned  Global  Hawk  aircraft 
taking  observations  over  the  large-scale  environment  from  an  altitude  that  provided  for 
profiles  of  the  entire  troposphere.  The  case  of  Typhoon  Sinlaku  is  used  as  operational 
typhoon  track  forecast  scenarios  were  highly  uncertain  and  research-based  aircraft 
observations  were  available  to  test  against  the  pseudo  Global  Hawk  observations.  Three 
numerical  experiments  are  conducted  such  that  initial  conditions  were  based  on 
conventional  observations,  conventional  and  inner-storm  aircraft  observations  and 
conventional  and  pseudo  Global  Hawk  observations.  Inclusion  of  pseudo  Global  Hawk 
observations  resulted  in  the  most  accurate  intensity  forecasts  but  track  forecasts  were  best 
using  observations  of  the  storm  inner  core. 
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I.  INTRODUCTION 


A,  HISTORY 

The  mission  of  U.  S.  Naval  Meteorology  and  Oceanography  Command  (METOC) 
has  long  been  centered  on  providing  warfighters  with  the  very  best  environmental 
information  they  need  to  make  the  most  informed  decision.  The  goal  is  to  never  be 
surprised  by  environmental  conditions.  One  of  the  most  severe  threats  to  naval  shore- 
based  and  afloat  resources  is  tropical  cyclones  (TCs),  which  occur  over  many  of  the 
ocean  basins  in  which  the  U.S.  Navy  operates.  As  such,  TC  reconnaissance,  analysis,  and 
forecasting  have  been  stated  by  USPACOM  as  essential  security  components  to  ensure 
safety  of  United  States  assets  and  National  Security  (Stephen,  2010). 

The  Joint  Typhoon  Weather  Center  (JTWC)  in  Pearl  Harbor,  Hawaii  is 
responsible  for  the  detection  and  prediction  of  the  path  and  strength  of  TCs  over  the 
western  north  Pacific  Ocean,  the  western  south  Pacific  Ocean,  and  the  north  and  south 
Indian  Oceans.  The  JTWC  mission  is  built  upon  the  need  to  protect  U.S.  assets,  maximize 
cost  efficiencies  by  balancing  risk  with  mission  goals;  and  supply  input  to  sortie  decision¬ 
making  processes.  To  achieve  these  goals,  the  most  accurate  estimates  of  current  and 
future  storm  positions  are  required.  In  addition  to  forecast  accuracy,  the  uncertainty  in 
forecast  positions  provides  key  information  under  which  a  decision  maker  must  utilize 
the  forecast.  As  forecast  uncertainty  decreases,  the  amount  of  geographic  area  that  may 
be  placed  into  unnecessary  warning  or  TC  conditions  of  readiness  (TCCOR)  will  be 
decreased. 

Much  earlier  than  the  advent  of  electronic  computing,  Lewis  Richardson  wrote  in 
1922  (Richardson  1965)  about  the  challenges  associated  with  numerically  integrating  the 
fundamental  equations  required  to  produce  a  very  rudimentary  time  evolution  of  500  hPa 
heights.  He  wrote,  “The  two  great  outstanding  difficulties  are  those  connected  with  the 
completeness  necessary  in  the  initial  observations  and  with  the  elaborateness  of  the 
subsequent  process  of  computing”  (Richardson,  1965).  The  basis  of  the  initial  conditions 
needed  to  initialize  a  time  integration  as  defined  by  Richardson  are  environmental 


1 


observations  that  did  not  become  routinely  available  until  much  later  than  the  seminal 
work  of  Richardson  (1965),  Lorenz  (1963)  defined  the  sensitivity  of  numerical 
simulations  to  small  variations  in  the  initial  conditions  upon  which  the  numerical 
integrations  were  begun.  Lorenz’s  discovery  of  large  variability  that  could  occur  in  the 
final  state  of  a  numerical  integration  when  initial  condition  variability  was  extremely 
small  lead  to  the  more  formal  concept  to  chaos  theory  (Gleick,  1987).  Lorenz’s  work  also 
lead  to  the  development  of  ensemble  forecasting  such  that  the  sensitivity  to  initial 
conditions  could  be  exploited  to  provide  an  estimate  of  the  uncertainty  or  predictability 
associated  with  a  given  numerical  forecast.  Eventually,  the  concept  contributed  to 
ensemble  data  assimilation  (Majumdar  and  et  ah,  2001),  which  can  be  used  to  define 
small,  random  perturbations  in  initial  conditions  based  on  estimates  of  observation 
accuracy  and  spatial  distribution.  In  this  thesis,  the  impact  of  data  distribution  on  initial 
conditions  and  subsequent  forecasts  is  defined  in  relation  to  the  problem  of  forecasting 
the  track  of  a  typhoon  over  the  western  North  Pacific. 

B,  MOTIVATION 

The  goal  of  this  thesis  is  to  build  on  the  program  titled  THe  Observing  System 
Research  and  Predictability  Experiment  (THORPEX)  and  more  specifically  the 
THORPEX-Pacific  Asian  Regional  Campaign  (T-PARC)  that  was  conducted  in  the  fall 
of  2008  (Elsberry  and  Harr  2008),  During  T-PARC,  more  than  1500  atmospheric 
soundings  were  captured  via  release  of  dropwindsondes  from  multiple  aircraft  platforms 
(Weissmann  et  ah,  2011).  While  the  aircraft  observations  were  obtained  to  provide 
observations  of  processes  (Raymond  and  Eopez-Carillo  2011,  Montgomery  et  al.  2010; 
Sanabia  2010)  important  to  TC  formation,  structure,  and  structure  change,  a  significant 
objective  of  T-PARC  was  to  examine  the  value  of  aircraft  observations  with  regard  to  the 
accuracy  of  TC  track  forecasts  produced  by  operational  numerical  forecast  models 
(Hamisch  and  Weissmann  2010;  Weissmann  et  al.  2011;  Chuo  et  al.  2012,  Weissmann  et 
al.  2012).  The  relative  value  of  observations  of  TC  location,  intensity,  structure,  and 
large-scale  environmental  conditions  were  obtained  over  regions  where  these 
observations  are  not  operationally  available. 
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During  T-PARC  on  11  September  2008,  a  three-plane  mission  (Figure  1)  was 
flown  into  Typhoon  (TY)  Sinlaku.  At  this  time,  there  was  considerable  uncertainty  in  the 
forecast  track.  For  this  aircraft  observing  period,  a  U.S.  Air  Force  WC-130J  aircraft 
conducted  a  standard  TC  reconnaissance  mission  into  the  storm  center.  Flight  level  for 
the  TC  reconnaissance  was  at  700  hPa  from  which  dropsondes  were  used  to  obtain 
profdes  of  winds,  temperature,  and  water  vapor.  The  stepped  frequency  microwave 
radiometer  (SFMR)  remotely  measured  surface  wind  speeds.  From  Taiwan,  the 
Dropsonde  Observations  for  Typhoon  Surveillance  near  the  Taiwan  Region  (DOTSTAR) 
aircraft  flew  near  200  hPa  and  deployed  dropwindsondes  in  the  immediate  vicinity  of  TY 
Sinlaku.  However,  the  DOTSTAR  aircraft  did  not  fly  over  the  center  of  the  storm. 
Finally,  the  Deutches  Zentrum  fur  Luft  (DLR)  FALCON  aircraft  flew  at  200  hPa  over  a 
region  poleward  of  Sinlaku.  The  purpose  of  the  FALCON  flight  was  to  gather 
observations  of  the  region  ahead  of  the  storm  and  measure  the  relative  strength  and 
positions  of  the  subtropical  ridge,  the  low-level  monsoon  trough,  and  an  upper-level 
midlatitude  trough. 

Hamisch  and  Weissmann  (2010),  examined  the  sensitivity  of  track  forecasts  of 
Sinlaku  (Figure  2)  produced  by  the  European  Center  for  Medium-Range  Weather 
Forecasts  (ECMWF)  numerical  forecast  model  to  variations  in  initial  conditions  based  on 
use  of  the  thee  aircraft  observation  sets  collectively  and  individually.  In  this  study,  they 
found  that  the  most  accurate  forecasts  were  produced  from  initial  conditions  in  which 
DOTSTAR  observations  from  high  altitude  in  the  near  storm  environment  were  used. 
These  forecasts  were  more  accurate  than  those  in  which  initial  conditions  were  produced 
by  assimilating  the  low-level  data  from  the  WC-130J  flight  in  the  typhoon  core  or  at  high 
altitude  in  a  remote  region  ahead  of  the  storm  as  observed  using  the  DER  Ealcon.  The 
above  example,  and  several  others  from  T-PARC  and  other  programs  provide  evidence  of 
the  value  of  in  situ  observations  of  a  TC  and  the  near  environment  of  the  TC  in  producing 
initial  conditions  that  lead  to  improve  forecasts  of  the  storm  position.  Eurthermore,  the 
results  from  T-PARC  point  to  the  variability  in  observation  value  based  on  location  of  the 
observation  relative  to  the  TC  and  to  the  information  with  respect  to  a  full  profile  of  the 
troposphere.  However,  the  only  ocean  basin  in  which  routine  aircraft  observations  are 
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taken  is  the  North  Atlantic.  Over  the  north  Atlantic,  the  U.S.  Air  Force  Weather 
Reconnaissance  Squadron  obtains  observations  in  the  TC  core  and  the  National  Oceanic 
and  Atmospheric  Administration  (NOAA)  Gulfstream  IV  obtains  observations  from  high 
altitudes  in  the  synoptic  environment  of  the  storm.  Given  the  sensitivities  to  observation 
location  and  altitude,  the  objective  of  this  thesis  is  to  examine  potential  improvements  to 
numerical  forecasts  of  TC  tracks  if  observations  were  available  from  high  altitude  in  the 
synoptic  environment  of  a  western  north  Pacific  typhoon,  where  routine  observations  are 
not  available. 


Figure  1.  Enhanced  infra-red  MTSAT  imagery  of  TY  Sinlaku  at  0915 
UTC  1 1  September  2008.  The  WC-130J  flight  track  is  defined 
by  the  black  line.  The  DOTSTAR  flight  track  is  defined  by  the 
red  line,  and  the  FALCON  flight  track  is  defined  by  the  yellow 
line,  (from  T-PARC  2008) 
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In  this  thesis,  an  experiment  is  eonducted  sueh  that  simulated  dropwindsondes  are 
defined  and  “deployed”  at  60,000  ft  as  would  be  eondueted  from  an  unmanned  aerial 
vehicle  (UAV)  such  as  the  Global  Hawk  (GH).  The  case  of  TY  Sinlaku  during  T-PARC 
is  chosen  for  the  experiment  because  there  is  a  unique  set  of  actual  aircraft  observations 
that  can  be  used  to  define  an  alternative  initial  condition  to  produce  a  type  of  benchmark 
forecast  track. 


Figure  2.  Track  forecasts  of  TY  Sinlaku  produced  by  ECMWF  global 
numerical  prediction  model  based  on  varying  initial  conditions. 
The  TC  positions  are  plotted  at  12  h  intervals.  Black  solid  dots 
are  best-track  locations.  Squares,  upward-pointing  triangles, 
downward-pointing  triangles,  and  circles  represent  forecast  of 
remote  observations,  vicinity  observations,  core  observations, 
and  all  observations,  respectively.  Cross  markers  represents 
forecast  of  no  observations  control  experiment,  (from  Harnisch 
and  Weissmann  2010) 


The  study  is  outlined  as  follows.  Chapter  II  details  the  methodology  of  the 
experiment  by  discussing  model  parameters,  experiment  design  and  execution.  Chapter 
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Ill  provides  a  background  by  describing  the  evolution  of  Typhoon  Sinlaku  in  September 
2008.  Chapter  IV  details  the  analysis  of  the  experimental  model  forecasts  by  comparing 
and  contrasting  the  various  forecasts  against  a  control  model  run  without  observation 
assimilation.  Chapter  V  provides  results  and  follow-on  study  opportunities. 


6 


II.  METHODOLOGY 


The  deterministic  and  ensemble  forecasts  of  TY  Sinlaku  were  performed  using 
the  Advanced  Weather  Research  and  Forecasting  (WRF-ARW)  non-hydrostatic 
mesoscale  model  (hereafter  WRF)  (Skamarock  et  al.  2008).  While  Sinlaku  existed  over 
the  western  North  Pacific  for  nearly  2  weeks  (Figure  3),  this  study  focuses  on  the  initial  5 
days  following  formation  of  the  storm.  During  this  period,  track  forecasts  were  highly 
variable  and  Sinlaku  underwent  significant  intensification  and  structure  changes. 

Two  nested  domains  (Figure  4)  with  horizontal  grid-spacing  of  27-  and  9-km  and 
42  vertical  levels  were  constructed  to  encompass  TY  Sinlaku.  The  near  environment 
microphysical  processes  were  handled  using  the  WRF  single  moment  6  class  (WSM6) 
microphysics  parameterization  scheme  (Hong  and  Lim  2006)  and  cumulus  convection 
was  accounted  for  using  the  Kain-Fritsch  scheme  (Fritsch  and  Kain  1993;  Kain  2004). 
The  Yonsei  University  (YSU)  boundary  layer  scheme  (Hong  et  al.  2006)  and  fifth 
generation  Penn  State  Mesoscale  Model  (MM5)  surface  layer  scheme  (Skamarock  et  al. 
2008)  were  used  to  represent  boundary  layer  and  surface  layer  processes,  respectively. 


Figure  3.  Sinlaku  track  comparison  where  the  black,  red,  green,  purple, 
and  blue  tracks  represent  JTWC  best  track,  ECMWF  YOTC 
analysis,  ECMWF  YOTC  forecast,  WRF  (YOTC)  analysis,  and 
WRF-(YOTC  forecast),  respectively.  (See  text  for  details) 
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Figure  4.  Domain  contiguration  used  tor  the  WKF  deterministic  torecasts 
and  the  DART-WRF  ensemble  data  assimilation  system.  The 
dOl  and  d02  labels  correspond  to  the  first  domain  (27-km  grid 
spacing)  and  second  domain  (9-km  grid  spacing),  respectively. 

Initially,  simulations  were  initialized  from  the  Global  Forecast  System  (GFS) 
reforecast  data  (Hamill  et  al.  2005),  however  the  initial  intensity  of  TY  Sinlaku  was 
unrealistically  weak  in  the  GFS  reforecast  data  (not  shown)  such  that  the  deterministic 
WRF  simulation  was  too  slow  to  develop  TY  Sinlaku  (not  shown).  In  addition,  the  GFS 
reforecast  data  were  limited  to  36  vertical  levels,  which  was  considered  too  coarse  to 
construct  pseudo-dropwindsonde  profiles  used  in  the  data  assimilation  aspect  of  this 
study.  Therefore,  the  initial  and  lateral  boundary  conditions  for  the  WRF  deterministic 
forecast  were  derived  from  the  ECMWF  operational  forecast  cataloged  in  the  Year  of 
Tropical  Convection  (YOTC,  Waliser  et  al.  2012)  data  archive.  Therefore,  these 
forecasts  initialized  at  1200  UTC  8  September  are  labeled  WRF-OPS.  It  is  important  to 
note  that  the  initial  conditions  for  the  WRF  forecasts  derived  from  the  ECMWF 
operational  forecast  represented  a  “best”  estimate  of  the  atmospheric  state  only  at  the 
initialization  time  of  the  global  forecast  (1200  UTC  8  September),  and  originated  from  a 
blend  of  the  previous  forecast  and  the  assimilation  of  observations  available  at  that  time. 
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A  second  WRF  forecast  was  generated  at  the  same  initial  time  (1200  UTC  8  September) 
but  boundary  conditions  were  defined  from  the  ECMWF  analysis  at  each  update  time. 
Therefore,  this  forecast  is  labeled  WRF-ANAL  and  has  information  not  available  in  real 
time.  The  WRF-ANAL  forecast  serves  as  a  type  of  upper-bound  in  accuracy  if  perfect 
boundary  information  would  be  available.  Figure  5  depicts  the  data  flow  described  above. 


Figure  5.  Data  flow  and  relationship  of  ECMWF  initial  and  boundary 

conditions  to  WRF  forecasts. 


To  determine  the  forecast  sensitivity  to  real  and  simulated  observations,  data 
assimilation  (DA)  was  conducted  using  the  Data  Assimilation  Research  Testbed  (DART) 
(Anderson  et  al.  2009)  ensemble  data  assimilation  system  in  conjunction  with  the  WRF 
model  (hereafter  referred  to  as  DART- WRF).  The  DART- WRF  ensemble  was  comprised 
of  48  members  that  used  an  identical  configuration  to  the  WRF  deterministic  forecasts 
described  previously. 

The  DA  cycle  (Figure  6),  which  is  more  generally  described  in  Anderson  et  al. 

(2009),  is  as  follows:  Observations  were  assimilated  for  the  period  1200  UTC  8 

September  -  0600  UTC  9  September  2008  using  a  3-h  assimilation  window  (±  90  minutes 
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centered  on  the  hour).  At  assimilation  time,  observations  within  the  ±  90  minute  window 
and  their  eorresponding  error  eharacteristics  were  eombined  with  the  ensemble  prior  state 
using  the  ensemble  adjustment  Kalman  filter  (EAKF)  (Anderson  2003)  to  generate  an 
increment  for  eaeh  ensemble  state.  Given  the  finite  ensemble  size,  whieh  leads  to 
sampling  error,  spatially  varying  state-space  inflation  (Anderson  2007)  was  also  used  to 
increase  spread  among  the  ensemble  members.  The  increment  was  then  applied  to  the 
ensemble  prior  state  to  generate  the  ensemble  posterior  state,  whieh  was  then  integrated 
forward  in  time  by  the  WRF  model  until  the  next  assimilation  time. 


3.  "Gross"  Obs.  Errors 


Figure  6.  Ensemble  data  assimilation  cycle  similar  to  that  used  for  the 
DART -WRF  simulations  (see  text  for  details).  The  model  is 
represented  by  the  green  arrows  and  the  state  vector  of  each 
ensemble  is  represented  by  the  blue  stars;  green  tiek  marks 
correspond  to  the  model  state  veetor  mapped  to  observation 
spaee  by  the  forward  operator  h;  red  tick  marks  represent 
observations,  and  red  eurves  correspond  to  the  observation 
probability.  The  blue  tiek  marks  represent  ensemble  updates 
provided  by  the  ensemble  filter;  blue  arrows  at  top  right  are  the 
increments  in  observation  space,  and  the  model  increments 
correspond  to  the  blue  arrows  near  bottom  center.  Sourees  of 
error  assoeiated  with  each  step  of  the  data  assimilation  eycle  are 
listed  in  the  figure  (from  Anderson  et  al.  2009). 
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Three  different  DA  experiments  were  eondueted  (Figure  7).  The  first  experiment, 
hereafter  WRF-CONV,  only  assimilated  eonventional  observations  that  were  available 
from  routine  data  sourees  (no  in  situ  observations).  The  second  DA  experiment,  hereafter 
WRF-C130,  assimilated  conventional  observations  as  well  as  dropwindsondes  deployed 
from  the  USAF  WC-130J  near  0300  UTC  9  September.  The  third  experiment,  hereafter 
WRF -PSEUDO,  assimilated  conventional  observations  and  pseudo  dropwindsondes 
generated  from  gridded  ECMWF  analysis  data  (Figure  7).  These  pseudo-dropwindsonde 
profiles  were  created  from  the  ECMWF  analysis  data  by  assuming  a  constant  fall  velocity 
of  12  m  s"'  from  a  height  of  18  km  (to  roughly  match  the  60,000  ft  cruising  altitude  of  the 
Global  Hawk)  while  allowing  for  horizontal  advection.  Trilinear  interpolation  was  used 
to  define  the  atmospheric  variables  at  the  pseudo-dropwindsonde  locations  based  on  the 
values  at  surrounding  ECMWF  grid  points.  Random  normal  noise  was  added  to  each 
profile  to  account  for  some  of  the  uncertainty  since  the  pseudo-dropwindsonde  profiles 
generated  from  the  gridded  model  analyses  were  much  ‘smoother’  than  actual 
dropwindsonde  profiles  (not  shown). 


DART-WRF 

WRF-OPS 

WRF-ANAL 


DART-WRF  ENSEMBLE  FORECASTS 


09/00 


10/00 


11/00 


12/00 


13/00 


time  (dd/hh) 


Figure  7.  Timeline  of  TY  Sinlaku  (2008)  forecast  experiments  conducted  in  this  study 
relative  to  available  in  situ  dropwindsondes  from  the  USAF  WC-130J  flight 
0133W  (red)  and  pseudo-dropwindsonde  observations  from  the  inner-core 
region  (green,  “PSEUDO- 1”)  and  surrounding  environment  (blue, 
“PSEUDO-2’’).  Numbers  in  parentheses  correspond  to  the  number  of 
dropwindsondes  for  each  observation  set.  The  red  and  green  lines  at  the  top 
of  the  figure  correspond  to  the  integration  times  of  the  WRF-ANAL  and 
WRF-OPS  deterministic  forecasts,  respectively.  The  blue  line  at  the  top  of 
the  figure  indicates  the  DART-WRF  ensemble  data  assimilation  phase  (1200 
UTC  8  September-0600  UTC  9  September),  and  the  purple  line  corresponds 
to  the  DART-WRF  ensemble  forecasts. 
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Since  ECMWF  analyses  were  only  available  every  6  h,  a  pseudo  flight-plan 
(Figure  8)  was  devised  based  on  an  assumed  cruising  speed  of  300  kt  to  fly  over  the  inner 
eore  (centered  at  0000  FITC  9  September)  and  through  the  outflow  of  the  storm  and  the 
midlatitude  jet  to  the  north  (eentered  at  0600  FITC  9  September).  Although  pseudo- 
dropwindsondes  were  generated  for  the  entire  flight  plan,  only  pseudo-dropwindsondes 
within  +/-90  minutes  of  0000  FITC  (red  dots  in  Figure  8a)  and  0600  FITC  (blue  dots  in 
Figure  8b)  9  September  were  assimilated.  Flnfortunately,  when  pseudo-dropwindsonde 
data  from  the  inner-eore  region  were  assimilated  during  the  0000  FITC  9  September 
assimilation  time,  over  half  of  the  ensemble  members  beeame  numerically  unstable  near 
the  observation  radius  of  influence  boundary.  This  may  have  been  a  result  of  the  data 
assimilation  conforming  too  heavily  to  the  pseudo  dropwindsonde  observations  near  the 
storm  center  such  that  there  was  a  large  disparity  between  the  region  impaeted  by  the 
pseudo-dropwindsondes  and  the  area  outside  the  influenee  of  these  observations.  As  a 
result,  only  pseudo-dropwindsonde  observations  from  0600  UTC  9  September  (Figure 
8b)  were  assimilated  for  the  WRF-PSFUDO  experiment. 

Observation  error  variances  for  all  observations  in  this  study  were  generated  from 
the  National  Center  for  Environmental  Prediction  (NCEP)  lookup-tables  based  on  the 
pressure  or  height  level  of  the  observation.  A  Gaspari-Cohn  eovarianee  cutoff  radius 
(Gaspari  and  Cohn  1999)  of  -630  km  was  used  for  all  observations,  which  limited  the 
observational  area  of  impact.  To  reduce  the  eomputational  expense,  dropwindsonde  data 
(pseudo  and  in  situ)  were  thinned  by  only  retaining  every  fifth  set  of  observations  in  the 
vertical. 

After  18  h  of  data  assimilation,  eaeh  ensemble  member  for  the  three  sets  of 
experiments  was  integrated  forward  in  time  from  its  state  at  0600  UTC  9  September  to 
0000  UTC  13  September  (for  an  additional  90  h).  To  analyze  the  impact  of  the  different 
observation  sets,  the  track,  intensity,  and  strueture  of  TY  Sinlaku  of  the  thee  ensemble 
means  were  examined  and  compared. 
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Figure  8.  Wind  speed  (m  s'^  shading  and  veetors)  and  geopotential  height 
(m,  eontours)  from  the  ECMWF  analysis  at  (A)  850  hPa  and  (B) 
200  hPa.  Panel  (A)  is  valid  at  0000  UTC  9  September  2008,  and 
panel  (B)  is  valid  at  0600  UTC  9  September  2008. The  blaek 
line  indicates  the  pseudo  flight  track  used  to  derive  the  pseudo- 
dropwindsonde  observations  (white  and  color-fdled  circles 
along  the  flight  track)  from  the  ECMWF  analysis.  The  red-filled 
dropwindsonde  locations  in  panel  (A)  indicate  the  inner-core 
pseudo-dropwindsondes  assimilated  at  0000  UTC  9  September, 
while  the  blue-fdled  dropwindsonde  locations  in  panel  (B) 
indicate  the  environmental  pseudo-dropwindsondes  assimilated 
at  0600  UTC  9  September.  The  starting,  middle,  and  ending 
times  (UTC)  are  indicated  for  each  pseudo-dropwindsonde  set. 
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III.  BACKGROUND 


A.  EVOLUTION  OF  SINLAKU 

Typhoon  Sinlaku  was  one  of  four  typhoons  to  be  observed  during  the  T-PARC 
experiment.  At  1800  UTC  8  September  2008,  Sinlaku  formed  just  west  of  Luzon  Island 
in  the  Philippines  (Figure  9)  The  T-PARC  focus  was  split  into  two  periods  during  the 
lifetime  of  Sinlaku.  A  series  of  aircraft  missions  were  conducted  beginning  on  9 
September  2008  (Chou  et  al.  201 1)  as  Sinlaku  moved  northward  toward  recurvature  over 
the  northern  portion  of  Taiwan.  Several  missions  were  then  conducted  following 
recurvature  when  Sinlaku  moved  into  the  midlatitude  and  underwent  extratropical 
transition  (Sanabia  2010). 


140'E 


Figure  9.  The  track  of  TY  Sinlaku  depicted  by  intensity  such  that  blue  is  a 


tropical  disturbance,  green  is  a  tropical  depression,  yellow  is  a 
tropical  storm,  red  is  a  typhoon,  and  magenta  is  an  extratropical 
cyclone.  The  numbers  in  each  circle  define  the  day  in 
September  2008.  (From  Kitamoto  National  Institute  of 
Informatics  (Nil)) 


At  1430  UTC  on  September  7  2008,  JTWC  issued  a  classification  of  “poor”  for 
an  area  of  interest  near  16. 3N  128. 5E  in  the  Philippine  Sea  due  to  persistent  convective 
activity  (JTWC,  2008).  Infrared  satellite  imagery  (Figure  10)  indicated  flaring  deep 
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convection  associated  with  a  developing  low-level  circulation  center  (LLCC).  At  this 
time,  the  indicated  wind  speed  from  the  advanced  scatterometer  (ASCAT)  ranged 
between  15  and  20  kt  (not  shown).  The  LLCC  was  developing  beneath  a  broad  upper- 
level  anticyclone.  Vertical  wind  shear  was  low,  while  a  moderate  level  of  upper  level 
diffluence  existed.  The  minimum  sea  level  pressure  (SLP)  was  estimated  to  be  near  1006 
hPa  with  winds  estimated  to  be  between  12  and  18  kt. 


Figure  10.  Enhanced  infrared  satellite  imagery  from  MTSAT2  at  1430 
UTC  7  September  2008.  The  black  arrow  points  to  an  area  of 
deep  convection  associated  with  the  pre-Sinlaku  disturbance. 

(From  T-PARC  Catalog) 

At  2200  UTC  7  September  2008,  JTWC  upgraded  the  classification  to  “fair”  for 
cyclone  development  over  the  next  24  hours  due  to  “favorable  upper  level  conditions 
despite  minimal  low  level  organization.”  Water  vapor  imagery  indicated  building  deep 
convection  over  a  weak  LLCC  (Figure  11).  The  strong  upper-level  outflow  on  the 
poleward  edge  of  the  convective  disturbance  was  enhanced  by  an  upper-level  cyclonic 
circulation  that  was  directly  to  the  east  of  the  pre-Sinlaku  disturbance.  At  this  time, 
JTWC  estimated  maximum  surface  level  winds  to  be  15  to  20  kt,  with  a  SLP  near  1005 
hPa. 
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Figure  11.  Water  vapor  imagery  MTSAT  at  2157  UTC  9  September,  2008 
showing  area  of  deep  convection.  (From  T-PARC  Catalog) 


At  0600  UTC  on  8  September  2008,  a  tropical  cyclone  formation  area  (TCFA) 
alert  was  issued  by  JTWC.  Deep  convection  continued  to  grow  near  a  consolidating 
LLCC.  An  upper  level  anticyclonic  circulation  along  with  low  vertical  wind  shear  and 
conditions  favorable  for  upper  level  diffluence  above  the  disturbance  area  were  all 
favorable  for  continued  development  of  the  convective  disturbance.  The  upper-level 
cyclonic  cell  to  the  east  continued  to  influence  the  outflow  on  the  poleward  side  of  the 
circulation.  Surface  winds  had  increased  to  near  23  kt  as  the  SLP  continued  to  fall  to 
1004  hPa.  JTWC  forecasters  upgraded  the  potential  for  a  significant  TC  to  occur  within 
the  following  24  hours  to  “good.” 

The  first  warning  was  issued  at  1500  UTC  8  September  for  tropical  depression 
15W.  Sustained  winds  had  reached  30  kt  with  gusts  of  40  kt.  The  consolidation  of  deep 
convection  at  the  storm  center  became  clear  in  infrared  imagery  taken  at  1430  UTC  8 
September  (Figure  12).  Deep  convective  cloud  bands  were  seen  wrapping  into  the  LLCC. 
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Figure  12.  Enhanced  infrared  satellite  imagery  from  MTS  AT  2  at  1430 
UTC  8  September  2008.  (From  T-PARC  Catalog) 

A  period  of  rapid  intensification  began  on  9  September  2008.  Over  the  previous 
period,  SEP  dropped  to  986  hPa,  and  surface  winds  increased  to  55-60  kt.  The  center  of 
vorticity  was  situated  beneath  a  well-defined  area  of  upper-level  divergence.  By  2032 
UTC  9  September,  a  WC-130J  flight  reported  a  SEP  of  938  hPa  and  a  maximum  surface 
wind  of  82  kt  as  measured  by  the  SFMR.  A  central  dense  overcast  (CDO)  was  becoming 
better  defined  with  intermittent  eye  wall  development  as  seen  in  visual  and  IR  imagery 
(not  shown).  The  outflow  jet  was  primarily  oriented  toward  the  northeast  quadrant  of  the 
storm.  The  development  of  Sinlaku  from  7-9  September  is  well  depicted  in  IR  imagery 
(Eigure  13). 

By  10  September,  a  well-defined  eye  was  visible  as  Sinlaku  reached  maximum 
intensity  with  a  SEP  of  927  hPa  and  maximum  winds  estimated  near  120  kt  (Eigure  14). 
In  satellite  imagery  there  is  evidence  of  dry  continental  air  to  the  north  and  west  of  the 
storm  (Eigure  13),  While  these  conditions  would  suggested  some  weakening  of  the  storm, 
the  storm  track  was  moving  slowly  in  a  northerly  direction  over  a  region  of  high  ocean 
heat  content  (OHC)  (Figure  15),  which  may  have  contributed  to  maintenance  of  the 
strong  intensity.  By  this  time,  the  previous  prediction  of  the  storm  to  recurve  had  been 
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reassessed  due  to  the  weakening  subtropical  steering  ridge  resulting  in  a  continued  drift 
toward  Taiwan. 


Figure  13.  Infrared  imagery  from  MTSAT  2  at  (a)  0230  UTC  7  September; 

(b)  0032  UTC  8  September  ;  (c)  0030  UTC  9  September;  and 
(d)  2030  UTC  9  September.  (From  T-PARC  Catalog) 
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Figure  14. 


Figure  15. 


Infrared  imagery  from  MTSAT2  at  1730  UTC  10  September 
(from  T-PARC  Catalog). 


Analyzed  OHC  (kJ  em'^)  valid  0000  UTC  10  September  2008. 
Cirele  represents  the  loeation  of  the  storm  (from  the  East  Asian 
Seas  Noweast  Foreeast  System  (EASNFS)). 
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An  eyewall  replacement  cycle  began  near  1750  UTC  10  September  (Figure  16a 
and  proceeded  over  the  next  17  h  to  1135  UTC  11  September  (Figures  16b-d).  The 
eyewall  replacement  cycle  coincided  with  some  weakening  in  storm  intensity.  Also  the 
storm  moved  over  a  region  of  decreased  OHC,  which  may  have  also  contributed  to  the 
decrease  in  storm  intensity  at  this  time  (Figure  17). 


Figure  16.  Eyewall  replacement  between  10  and  1 1  Sept  08  (after 

TPARC_2008  weather  summary) 

Over  the  next  24  hours,  Sinlaku  tracked  to  the  northwest  and  continued  to 
weaken.  The  weakening  was  likely  due  to  increased  vertical  wind  shear  (not  shown)  and 
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movement  over  a  region  of  reduced  OHC  (Figure  17).  At  1750  UTC  12  September, 
aircraft  missions  deploying  dropsondes  reported  maximum  surface  winds  of  98  kt  and 
SLP  of954hPa. 


Figure  17.  Analyzed  OHC  (kJ  cm'^)  valid  1800  UTC  12  September  2008 
(from  T-PARC  Catalog  and  Naval  Research  Lab  East  Asian 
Seas  Nowcast  /  Forecast  System  (NRL  EASNFS). 

As  Sinlaku  approached  Taiwan  during  12-13  Sept,  the  eyewall  contracted,  which 
is  clearly  identified  in  microwave  imagery  (Eigure  18).  At  this  time,  the  satellite-based 
intensity  estimate  was  102  kts  vertical  wind  shear  increased  (not  shown)  and  Sinlaku 
made  landfall  on  13  September.  Reported  observations  from  the  Yonaginijima  Island 
(eye  did  not  pass  directly  over  the  island)  reported  maximum  winds  near  92  kt  and  SLP 
of  962  hPa  (Eigure  19). 

Over  the  ensuing  days,  Sinlaku  continued  to  weaken  and  moved  slowly  to  the  NW 
before  recurvature  on  14  September.  Comparison  of  IR  imagery  comparisons  from  13-14 
Sept  showed  considerable  decreased  intensity. 
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Figure  18.  Combined  85  GHz  microwave  and  infrared  imagery  at  (a)  2157 
UTC  12  September  and  (b)  0941  UTC  13  September  (from 
Naval  Research). 


13  September  2008  Yonagunijima 
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Figure  19.  Observations  from  Yonaginijima  Airport  13  Sep  08  (from 

TPARC_2008  weather  summary) 
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After  passing  by  Taiwan  and  recurving,  Sinlaku  nearly  became  stationary  for  a 
period  of  48  hours.  By  18  September,  another  period  of  rapid  intensification  occurred  as 
winds  increased  by  20  knots  over  an  18-h  period  and  Sinlaku  regained  typhoon  intensity. 
Finally  on  20  September,  2008,  a  final  warning  was  issued  on  Sinlaku  as  it  moved  east  of 
Tokyo,  Japan.  During  the  lifecycle  of  Sinlaku,  47  warnings  were  issued  by  JTWC,  and 
the  max  sustained  winds  topped  out  at  125  knots. 

B,  SYNOPTIC  ANALYSIS  DURING  SINLAKU 

The  synoptic  environment  of  the  pre-recurvature  period  of  Sinlaku  is  defined 
using  the  National  Centers  for  Environmental  Prediction  (NCEP)  Climate  Eorecast 
System  Reanalysis  (CESR)  (Saha  et  al.  2010).  The  NCEP  CESR  has  38  km  global 
resolution  and  fields  available  every  6  h.  The  primary  synoptic  factors  relevant  to  the 
movement  of  Sinlaku  are  identified  for  eventual  comparison  among  the  experimental 
forecast  fields  initialized  using  varying  data  sources. 

The  synoptic  analysis  identifies  the  basic  large-scale  features  over  the  western 
North  Pacific  that  influenced  the  motion  of  Sinlaku  during  this  period  of  high  forecast 
uncertainty.  The  role  of  the  low-level  monsoon  trough  and  subtropical  ridge  is  examined. 
The  role  of  a  midlatitude  trough  is  also  examined  in  relation  to  upper-level  outflow 
patterns  and  the  influence  on  recurvature. 

1.  Sea  Level  Pressure 

During  the  transition  season  in  the  western  north  Pacific  the  low-level  atmosphere 

can  be  generally  characterized  as  a  near  equatorial  monsoon  trough  that  extends  from  the 

west  to  the  Philippines  with  temporary  eastward  extension  over  the  Philippine  Sea. 

Convection  often  forms  in  the  monsoon  trough  and  contributes  large-scale  upper  level 

mass  divergence.  At  the  eastern  edge  of  the  monsoon  trough,  easterly  trade  winds 

crossing  the  Pacific  and  meet  the  westerly  winds  associated  with  the  equatorial  portion  of 

the  monsoon  trough.  This  defines  a  region  of  low-level  confluence  and  convergence  that 

contributes  to  the  forcing  of  deep  convection  that  may  propagate  westward  along  the 

monsoon  though  axis.  Just  poleward  of  the  monsoon  trough,  a  strong  subtropical  ridge 

covers  most  of  the  subtropical  western  north  Pacific.  The  strong  easterly  flow  between 
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the  poleward  portion  of  the  monsoon  trough  and  the  equatorward  side  of  the  subtropieal 
ridge  provide  a  general  steering  flow  in  whieh  most  TCs  that  form  in  the  trough  move 
northwestward.  If  the  subtropieal  ridge  is  weakened  by  the  influence  of  a  midlatitude 
trough,  the  general  easterly  steering  flow  weakens  and  a  TC  may  become  influence  by 
the  westerly  midlatitude  flow  and  recurve  into  the  midlatitudes. 

Sinlaku  formed  within  the  low-level  trough  associated  with  the  western  North 
Pacific  summer  monsoon  (WNPSM).  Throughout  the  2008  summer,  WNPSM  was  weak 
as  trade  winds  dominated  the  region  of  the  Philippine  Sea.  Enhanced  areas  of  active 
monsoonal  flow  during  the  onset  of  Sinlaku  in  early  September  2008  had  the  effect  of 
increased  convective  activity,  clouds  and  precipitation.  On  1200  UTC  9  September,  the 
low  pressure  center  for  Sinlaku  was  situated  within  the  WNPSM  just  east  of  Luzon, 
Philippines  (Figure  20).  A  weak  low  pressure  system  to  the  northeast  of  Sinlaku  was  a 
low-level  reflection  of  an  upper-level  trough.  The  combination  of  the  Sinlaku  disturbance 
and  the  low  pressure  to  the  northeast  contributed  to  a  weakened  subtropical  ridge  and  the 
slow  northward  movement  of  Sinlaku. 


Figure  20.  Analyzed  mean  sea  level  pressure  (contours  in  2  hPa 
increments)  for  1200  FITC  9  September  2008 


25 


2. 


850  hPa 


The  subtropical  ridge  was  pushed  farther  poleward  than  normal  due  to  the 
orientation  of  the  monsoon  trough  to  the  northeast  from  Sinlaku  to  the  low  pressure  over 
the  subtropics  northeast  of  Sinlaku  (Figure  21).  These  low  pressure  systems  forced  the 
western  extension  of  the  subtropical  ridge  to  be  over  Japan. 


Figure  21 .  Analyzed  850  hPa  winds  (m  s’')  and  geopotential  height 

(contours,  in  30  m  intervals)  for  1200  UTC  9  September  2008 

The  forcing  of  the  subtropical  ridge  to  the  north,  reduced  the  easterly  steering 
flow  over  the  environment  of  Sinlaku  and  resulted  in  the  slow  northward  movement  of 
the  TC  during  this  time  period. 
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3. 


500  hPa 


An  upper-level  midlatitude  trough  was  moving  off  of  eastern  China  and  began  to 
erode  the  subtropieal  (figure  22).  As  Sinlaku  became  more  mature  and  moved  northwest, 
its  speed  of  advance  was  too  slow  and  Sinlaku  was  not  caught  up  by  steering  to  the  east 
due  to  the  long  wave  trough  to  the  north.  This  would  have  presumably  caused  Sinlaku  to 
recurve  sooner.  Consequently,  the  turn  of  Sinlaku  to  the  north  and  eventual  northeast  was 
delayed  until  after  landfall  over  northern  Taiwan. 

4.  200  hPa 

At  200  hPa,  the  outflow  from  Sinlaku  was  turned  sharply  to  the  east  and 
equatorward.  The  turning  of  the  outflow  was  due  to  the  influence  of  the  upstream 
midlatitude  trough.  The  resulting  strong  confluent  westerly  flow  at  200  hPa  was 
immediately  poleward  of  the  storm  (Figure  22). 


Figure  22.  Analyzed  mean  sea  level  pressure  at  1200  UTC  9  September 
2008.  Dashed  and  solid  lines  represent  500  hPa  and  200  hPa 
geopotential  heights  (m),  respectively.  Wind  barbs  are  at  200 
hPa  with  small  barbs  representing  5  m  s'^ 
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IV.  OBSERVATION  SENSITIVITY  EXPERIMENTS 


A,  DETERMINISTIC  WRF  MODEL  (WRF-OPS)  SIMULATIONS 

The  initial  and  lateral  boundary  conditions  for  the  WRF  deterministic  forecast 
WRF-OPS  and  the  WRF  ensemble  forecasts  were  derived  from  the  ECMWF  operational 
forecast  initialized  at  1200  UTC  8  September  2008.  To  simulate  real-time  operations,  the 
experimental  model  forecasts  utilized  the  ECMWF  forecast  (labeled  ECMWF -YOTC) 
for  1200  UTC  8  September  to  define  the  WRF-OPS  single,  deterministic  model  forecast. 
No  data  assimilation  was  conducted  in  the  production  of  the  WRF-OPS  forecast.  This 
allowed  for  comparison  with  the  experimental  runs  to  follow:  WRF-CONV,  WRF-C130, 
and  WRF-PSEUDO,  which  all  included  data  assimilation. 

The  forecast  track  and  intensity  forecast  of  WRF-OPS  were  compared  with  the 
JTWC  best  track  data  (Figure  23).  In  the  WRF-OPS  forecasts,  Sinlaku  made  landfall  on 
the  southern  end  of  Taiwan  whereas  the  actual  storm  crossed  the  very  northern  portion  of 
Taiwan.  In  contrast,  the  WRF-ANAL  forecast  was  more  in  agreement  with  JTWC  best 
track  (Figure  23),  as  expected.  The  striking  departure  of  WRF-OPS  from  the  best  track 
and  WRF-ANAE  suggests  a  cumulative  positive  effect  of  observations  on  TC  forecasts. 
The  forecast  intensity  from  WRF-ANAL  and  WRF-OPS  (Figure  24)  matched  the  JTWC 
best  track  intensity  except  near  1 1  and  12  September.  At  this  time  in  WRF-OPS  forecasts, 
Sinlaku  was  overland  in  southern  Taiwan  so  the  mean  sea-level  pressure  (MSLP) 
increased  rapidly  (Figure  24).  Both  WRF  intensity  forecasts  are  much  more  accurate  than 
the  operational  ECMWF  forecasts  or  the  analyzed  ECMWF  intensities. 


29 


Figure  23.  Sinlaku  track  comparison  where  the  black,  red,  green,  purple, 
and  blue  tracks  represent  JTWC  best  track,  ECMWF  YOTC 
analysis,  ECMWF  YOTC  forecast,  WRF  (YOTC)  analysis,  and 
WRF-(YOTC  forecast),  respectively. 


During  this  period,  JTWC  forecast  Sinlaku  to  recurve  to  the  east  much  sooner 
than  it  actually  did.  The  JTWC  warning  number  four  initialized  at  0600  UTC  on  9 
September  forecasted  a  304  nm  closest  point  of  approach  (CPA)  for  Taipei  (Figure  25). 
Just  48  hours  later,  CPA  to  Taipei  was  reduced  to  92  nm  (Figure  26).  By  14  September 
the  eye  of  Sinlaku  was  passing  over  Taipei  (not  shown).  The  uncertainty  of  the  forecast 
was  well  represented  by  the  large  area  of  possible  track  locations  defined  in  relation  to 
each  forecast  (shaded  region  in  Figures  25  and  26). 
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Figure  24.  Sinlaku  intensity  eomparison  where  the  black,  red,  green, 
purple,  and  blue  lines  represent  JTWC  best  track,  ECMWF 
YOTC  analysis,  ECMWF  YOTC  forecast,  WRF-ANAE  and 
WRF-OPS,  respectively. 
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Figure  25.  JTWC  TC  warning  number  4  from  0600  UTC  9  September 

2008  (from  JTWC). 
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Figure  26.  JTWC  TC  warning  number  12  from  0600  UTC  1 1  September 

2008  (from  JTWC). 

The  WRF-OPS  forecast  verifying  at  1200  UTC  8  September  (Figure  27)  had  a 
fairly  weak  ridge  at  850  hPa  north  and  east  of  Sinlaku.  This  is  defined  by  the  red  1530  m 
contour  in  Figure  27.  A  well-established  monsoon  trough  extended  from  the  storm  center 
to  the  southwest  (yellow  contour  in  Figure  28).  A  stronger  ridge  existed  at  500  hPa, 
which  is  defined  by  the  red  contour  in  Figure  28.  The  most  prominent  upper  level  feature 
at  200  hPa  (Figure  29)  is  the  outflow  that  extends  north  of  the  storm.  The  confluence  of 
the  outflow  and  the  eastward-moving  trough  leads  to  strong  westerly  flow  north  of 
Sinlaku  (Figure  29). 

In  the  12  h  forecast  (Figure  30)  verifying  at  0000  UTC  9  September,  the  850  hPa 
ridge  strengthened  and  pushed  farther  south  toward  Sinlaku.  The  storm  intensified  and 
the  monsoon  flow  that  encompassed  the  TC  had  intensified.  As  the  ridge  strengthened, 
the  storm  moved  more  northwestward  as  defined  by  the  red  motion  vector  in  Figure  30. 
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Figure  27.  Analyzed  winds  (veetors)  and  speed  (m  s'\  shaded)  at  850  hPa 
and  850  hPa  geopotential  heights  (m,  eontours  at  20  m 
intervals).  The  red  veetor  represents  general  storm  motion.  The 
red  and  yellow  eontours  represent  1530  m  and  1510  m 
geopotential  heights,,  respeetively,  at  1200  UTC  8  September. 
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Figure  28.  Analyzed  winds  (veetors)  and  speed  (m  s'\  shaded)  at  500  hPa 
and  500  hPa  geopotential  heights  (m,  eontours  at  20  m  intervals) 
The  red  and  yellow  eontours  represent  5880  m  and  5860  m 
geopotential  heights,,  respeetively,  at  1200  UTC  8  September. 
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Figure  29.  Analyzed  winds  (vectors)  and  speed  (m  s'*,  shaded)  at  200  hPa 
and  200  hPa  geopotential  heights  (m,  contours  at  20  m  intervals) 
The  red  and  yellow  contours  represent  12,460  m  and  12,440  m 
geopotential  heights,  respectively,  at  1200  UTC  8  September. 
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Figure  30.  WRF-OPS  12-h  forecast  winds  (vectors)  and  speed  (m  s'^ 
shaded)  at  850  hPa  and  850  hPa  geopotential  heights  (m, 
contours  at  20  m  intervals).  The  red  vector  represents  general 
storm  motion.  The  red  and  yellow  contours  represent  1520m 
and  1500  m  geopotential  heights,  respectively,  at  0000  UTC  9 

September. 
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In  the  18  h  forecast  that  verities  at  0600  UTC  on  09  September  the  surrounding 
environment  of  Sinlaku  continued  to  change.  The  850  hPa  ridge  (red  1520  m  contour 
Figure  31)  moved  poleward.  At  500  hPa  the  subtropical  ridge  defined  by  the  red  5880  m 
contour  in  Figure  32  weakened  north  of  Sinlaku  was  positioned  just  north  of  the  storm. 
The  ridge  was  no  longer  continuous  north  of  the  storm,  which  had  begun  to  turn  toward 
the  north.  At  200  hPa  (Figure  33),  the  anticyclone  over  the  TC  expanded.  To  the  north  of 
the  TC  the  midlatitude  trough  developed  a  closed  circulation  just  southwest  of  Japan.  The 
combination  of  the  outflow  and  midlatitude  trough  contributed  to  a  strengthening  of  the 
westerly  flow  (Figure  33)  north  of  the  storm. 
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Figure  3 1 .  WRF-OPS  12-h  forecast  winds  (vectors),  and  speed  (m  s’', 
shaded)  at  850  hPa  and  850  hPa  geopotential  heights  (m, 
contours  at  20  m  intervals)  The  red  and  yellow  contours 
represent  1520m  and  1500  m  geopotential  heights,  respectively, 
at  0600  UTC  9  September. 
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Figure  32.  WRF-OPS  12-h  forecast  winds  (vectors)  and  speed  (m  s'^  shaded)  at 
500  hPa  and  500  hPa  geopotential  heights  (m,  contours  at  20  m 
intervals)  The  red  and  yellow  contours  represent  5880  m  and  5860  m 
geopotential  heights,  respectively,  at  0600  UTC  9  September. 
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Figure  33.  WFF-OPS  12-h  forecast  winds  (vectors)  and  speed  (m  s'\  shaded)  at 
200  hPa  and  200  hPa  geopotential  heights  (m,  contours  at  20  m 
intervals)  The  red  and  yellow  contours  represent  12,460  m  and  12,440 
m  geopotential  heights,  respectively,  at  0600  UTC  9  September. 
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At  1200  UTC  on  9  September,  the  large-scale  environment  continued  to  evolve  as 
Sinlaku  tracked  very  slowly  to  the  northwest.  At  850  hPa  (Figure  34),  the  western 
portion  of  the  ridge  to  the  north  of  the  TC  weakened.  The  weak  mid-level  ridge  (red  5880 
m  contour  in  Figure  35)  that  previously  stretched  north  of  Sinlaku  expanded  over  a  larger 
area  to  the  east  of  the  storm  (Figure  35).  Flowever,  the  ridge  to  the  north  of  Sinlaku  was 
still  weak  as  the  upper-level  cyclone  circulation  evolved  at  200  hPa  at  0600  UTC  (Figure 
35)  now  extended  down  to  500  hPa.  At  200  hPa  level  (Figure  36),  the  ridge  over  the 
storm  intensified.  The  outflow  from  Sinlaku  continued  to  merge  with  the  digging  and 
eastward-moving  trough,  which  resulted  in  increased  westerly  winds  to  the  north  of  the 
storm. 
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Figure  34.  WRF-OPS  24-h  forecast  winds  (vectors)  and  speed  (m  s'\ 
shaded)  at  850  hPa  and  850  hPa  geopotential  heights  (m, 
contours  at  20  m  intervals).  The  red  vector  indicates  general 
storm  motion.  The  red  and  yellow  contours  represent  1520  m 
and  1500  m  geopotential  heights,  respectively,  at  1200  UTC  9 

September. 
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Figure  35.  WRF-OPS  24-h  forecast  winds  (vectors)  and  speed  (m  s'^  shaded)  at  500 
hPa  and  500  hPa  geopotential  heights  (m,  contours  at  20  m  intervals).  The 
red  contour  represents  5880  geopotential  heights  at  1200  UTC  9  September 
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Figure  36.  WRF-OPS  24-h  forecast  winds  (vectors)  and  speed  (m  s'^  shaded)  at 
200  hPa  and  200  hPa  geopotential  heights  (m,  contours  at  20  m 
intervals).  The  red  and  yellow  contours  represent  12,480  m  and  12460  m 
geopotential  heights,  respectively,  at  1200  UTC  9  September. 
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In  the  36  h  forecast  that  verified  at  0000  UTC  10  September  (Figure  37),  the  low- 
level  ridge  (red  1530  contour  in  Figure  37)  continued  to  weaken.  The  most  dramatic 
change  occurred  at  mid-levels  (Figure  38)  as  the  subtropical  ridge  north  of  the  storm  was 
replaced  by  the  deepening  midlatitude  trough.  However,  the  ridge  did  extend  southward 
to  be  northeast  of  Sinlaku  (Figure  38).  The  orientation  was  likely  responsible  for  steering 
the  storm  more  westward  toward  southern  Taiwan.  As  the  storm  matured,  the  200  hPa 
heights  directly  above  the  storm  continued  to  increase  (Figure  39).  The  outflow  from 
Sinlaku  continued  to  increase  and  wrap  around  to  the  east  and  south  (Figure  39).  The  mid 
and  lower  level  ridges  continued  to  move  east. 
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Figure  37.  WRF-OPS  36-h  forecast  winds  (vectors)  and  speed  (m  s'^ 
shaded)  at  850  hPa  and  850  hPa  geopotential  heights  (m, 
contours  at  20  m  intervals).  The  red  vector  indicates  general 
storm  motion.  The  red  and  yellow  contours  represent  1520  m 
and  1500  m  geopotential  heights,  respectively,  at  0000  UTC  10 

September. 
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Figure  38.  WRF-OPS  36-h  forecast  winds  (vectors)  and  speed  (m  s-1,  shaded)  at 
500  hPa  and  500  hPa  geopotential  heights  (m,  contours  at  20  m 
intervals).  The  red  and  yellow  contours  represent  5885  m  and  5865  m 
geopotential  heights,  respectively,  at  0000  UTC  10  September. 


Figure  39. 


WRF-OPS  36-h  analyzed  winds  (vectors),  speed  (m  s'\  shaded)  at  200  hPa 
and  200  hPa  geopotential  heights  (m,  contours  at  20  m  intervals).  The 
yellow  contour  represents  12,460  m  geopotential  heights  at  0000  UTC  10 

September. 
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Over  the  next  36  hours,  the  72  h  forecast  that  verified  at  1200  UTC  1 1  September 
indicated  that  landfall  on  the  southern  end  of  Taiwan  was  imminent.  Under  the  influence 
of  the  strong  midlatitude  trough,  lower-level  (Figure  40)  and  mid-level  (Figure  41)  ridge 
retreated  poleward  over  Japan.  The  upper  level  trough  (Figure  42)  was  moving  east  of 
Japan.  As  the  trough  moved  eastward,  its  influence  to  Sinlaku  weakened  as  Sinlaku 
continued  to  move  west-northwest  toward  southern  Taiwan. 
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Figure  40.  WRF-OPS  72-h  forecast  winds  (vectors)  and  speed  (m  s'^ 
shaded)  at  850  hPa  and  850  hPa  geopotential  heights  (m, 
contours  at  20  m  intervals).  The  red  vector  indicates  general 
storm  motion.  The  red  and  yellow  contours  represent  1500  m 
and  1480  m  geopotential  heights,  respectively,  at  1200  UTC  1 1 

September. 
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Figure  41 .  WRF-OPS  72-h  forecast  winds  (vectors)  and  speed  (m  s’',  shaded) 
and  500  hPa  geopotential  heights  (m,  contours  at  20  m  intervals).  The 
red  and  yellow  contours  represent  5880  m  and  5860  m  geopotential 
heights,  respectively,  at  1200  UTC  11  September. 
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Figure  42.  WRF-OPS  72-h  forecast  winds  (vectors)  and  speed  (m  s’',  shaded)  at 
200  hPa  and  200  hPa  geopotential  heights  (m,  contours  at  20  m 
intervals).  The  red  and  yellow  contours  represent  12,480  m  and  12,460 
m  geopotential  heights,  respectively,  at  1200  UTC  11  September. 
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As  indicated  in  Figure  3,  the  WRF-OPS  forecast  track  of  Sinlaku  was  too  far 
south  than  the  aetual  traek.  Based  on  the  analysis  of  the  WRF-OPS  foreeast  fields,  the 
southward  traek  was  due  to  the  subtropieal  ridge  remaining  a  strong  infiuenee  to  the 
northeast  of  the  storm.  Also,  the  midlatitude  trough  moved  rapidly  eastward  sueh  that  it 
did  not  infiuenee  Sinlaku. 

B,  SINLAKU  WRF-DA  EXPERIMENTS 

In  this  study,  three  primary  parameters  (storm  intensity,  steering  flow,  and  initial 
position  of  the  storm)  were  examined  to  determine  the  impaets  of  observation  and  type 
were  on  the  model  foreeasts  of  Sinlaku.  During  Sinlaku,  the  operational  foreeast 
uncertainty  assoeiated  with  storm  traek  and  intensity  was  high.  The  time  and  loeation  of 
reeurvature  to  the  north  then  east  was  not  well  foreeasted  by  the  numerieal  weather 
models.  It  is  widely  held  that  a  better  initial  analysis  of  the  large-scale  eonditions  around 
the  storm  often  lead  to  a  more  aecurate  foreeast.  In  an  attempt  to  gain  more  confidenee  in 
the  forecast,  this  study  examined  potential  impacts  of  strategic  observations.  To  this  end, 
thee  WRF  experimental  model  runs  were  produeed  and  examined  to  see  what  impaets 
were  made. 

Various  sets  of  observations  were  assimilated  to  define  the  WRF  model 
experiments.  As  diseussed  in  Chapter  II,  the  first  variant  (WRF-CONV)  assimilated  only 
eonventional  observations  into  the  model.  The  seeond  variant  (WRF-C130)  assimilated 
both  conventional  observations  as  well  as  the  aetual  WC-130J  dropwindsonde  data 
eolleeted  during  T-PARC  and  Typhoon  Sinlaku.  The  final  variation  of  the  WRF  model 
run  (WRF-PSEUDO)  assimilated  eonventional  observations  and  pseudo  dropwindsonde 
data  derived  from  ECMWF  analysis  of  Sinlaku. 

1.  Initial  Position 

Many  faetors  ean  be  attributed  to  the  differenees  among  numerieal  foreeasts.  Eor 
example,  storm  strueture  and  intensity  ean  play  a  role  as  well  as  synoptie  weather 
patterns.  In  this  study,  it  is  also  important  to  note  that  eaeh  WRE  simulation  eonducted 
was  by  design  assimilated  with  different  sets  of  observation.  The  result  not  only 
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produced  different  ensemble  consensus,  but  also  established  different  initial  positions  of 
the  storm  center 

The  experimental  WRF  ensemble  model  runs  (Figure  43)  indicate  various  degrees 
of  uneertainty  and  accuracy.  The  degree  of  uncertainty,  whieh  is  based  on  the  spread 
among  ensemble  members,  between  WRF-CONV  and  WRF-PSEUDO  is  very  similar. 
However,  there  are  subtle  differences.  The  WRF-CONV  is  biased  toward  a  more 
southerly  landfall  position  while  the  WRF-PSEUDO  traek  ensemble  average  has  a  more 
northerly  track.  The  traek  of  WRE-CONV  is  closest  to  the  WRE-OPS,  whieh  indieates 
that  the  assimilation  of  eonventional  observations  is  likely  similar  to  the  use  of  the 
ECMWE  operational  foreeasts.  This  is  eonsistent  as  the  ECMWE  foreeast  were 
initialized  with  the  conventional  observations.  Both  WRE-CONV  and  WRE-PSEUDO 
contain  more  track  uncertainty  than  does  WRE-C130.  It  is  hypothesized  that  the  WRE- 
C130  had  the  added  benefit  of  21  in  situ  observations  of  the  storm  center  as  well  as 
conventional  observation  assimilation.  The  in  situ  observations  of  the  storm  center  meant 
the  aecuracy  of  the  storm  plaeement  during  the  time  of  the  model  initialization  was  much 
greater  than  that  of  WRF-CONV  and  WRF-PSEUDO.  By  having  a  better  sense  of  the 
eurrent  storm  structure,  the  WRE-C130  foreeast  produeed  less  uncertainty  and  a  more 
accurate  forecast  than  did  the  model  with  conventional  observations  alone. 


44 


WRF-CONV  WRF-C130  WRF-PSEUDO 


Figure  43.  Ensemble  WRF  foreeasts  from  the  three  WRF-DA  experiments. 

Grey  lines  are  ensemble  member  forecast  tracks  for  TY  Sinlaku. 

Orange  line  is  the  ensemble  consensus.  Black  line  is  JTWC 
best  track.  Red  line  is  ECMWF  YOTC  analysis.  Green  line  is 
ECMWF  YOTC  forecast.  Purple  line  is  WRF-ANAL.  Blue 
line  is  WRF-OPS  forecast  (no  additional  observation 
assimilation).  Dots  represent  potential  storm  position  for  the 
given  model  in  24  h  intervals. 

Although  WRF-PSEEfDO  did  not  lead  to  a  more  accurate  forecast  over  WRF- 
C130,  it  was  improved  over  WRF-CONV.  In  comparison,  WRF-PSEUDO  included 
conventional  observations  and  ECMWF  derived  pseudo  dropwindsonde  observations  of 
the  upper-level  outflow  jet  region  of  the  storm.  With  only  15  additional  observations 
along  the  outflow  jet,  the  resulting  forecast  shifted  what  was  a  southerly  track  of  WRF- 
CONV  to  a  more  northerly  track  of  WRF-PSEUDO,  thus  bringing  the  forecast  in  greater 
alignment  with  the  actual  track. 

Because  the  WC-130J  observations  were  obtained  at  low-levels  in  the  storm  core, 
it  is  expected  that  assimilation  of  those  data  would  yield  a  more  accurate  storm  structure 
and  initial  position.  Comparison  of  storm  position  between  the  WRF-C130  and  WFD- 
CONV  (Figure  44)  defines  a  vorticity  dipole,  which  is  indicative  of  differences  in 
position.  The  dipole  pattern  suggests  that  the  storm  was  moving  too  quickly  to  the  north 
and  east  in  the  WRF-CONV  forecast.  Since  the  experimental  models  started  the  storm  in 
different  locations  and  the  models  were  built  on  different  observation  sets,  the  resulting 
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ensemble  predictions  were  also  different.  The  degree  to  which  initial  positional 
difference  played  in  each  ensemble  consensus  track  is  difficult  to  state,  however  it  can  be 
speculated  that  the  initial  position  of  the  storm  is  a  large  contributor  to  the  resulting 
accuracy  of  the  forecast. 


2008-09-09-1500  average  relative  vortidly  850  hPa  (lO"®  s'') 

SINLAKU_2  (48)  -  SINLAKU_1  (48)  average  geopotential  height  (m) 


Figure  44.  Vorticity  difference  WRF-C130  minus  WRF-CONV  for  the  18 
h  forecast  verified  at  1500  UTC  9  September.  Shading 
represents  relative  vorticity  differences  at  850  hPa  [10'^  s'*]. 
Contours  represent  geopotential  height  differences  [m]. 


It  is  possible  that  if  storm  centered  in  situ  observations  had  been  successfully 
assimilated  into  WRF-PSEUDO,  the  resulting  ensemble  forecast  would  have  been  more 
accurate  thank  the  WRF-C130. 

2,  Mean  Wind  Steering  Flow 

The  steering  flow  defined  in  each  WRF  model  forecast  was  analyzed  to  help 
explain  the  forecast  track  differences.  Wind  difference  fields  were  also  examined 
between  WRF  forecasts  to  ascertain  causes  for  the  steering  flow  changes  from  forecast  to 
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forecast.  Holland  (1984)  defined  steering  flow  as  a  weighted  average  wind  between  300 
and  800  hPa  pressure  heights  extended  over  5  to  7  degrees  latitude  from  the  center  of 
rotation.  Fiorino  and  Elsberry  (1989)  went  further  by  advocating  a  larger  annulus  size  of 
300-1000  km  from  the  storm  center  for  averaging  steering  flow.  The  steering  flow 
calculations  in  this  study  used  various  combinations  of  pressure  levels  and  annulus  sizes 
to  generate  several  steering  flow  definitions.  Among  them  were  the  parameters  advocated 
by  Holland  (1984)  and  Fiorino  and  Elsberry  (1989).  All  plots  were  examined  carefully  to 
determine  which  levels  and  annulus  size  bore  the  greatest  resemblance  to  the  average 
storm  motion  plot  of  each  WRE  run  (Eigure  45).  Then  the  relevant  atmospheric  factors 
were  considered  in  each  WRE  solution. 
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Eigure  45.  Average  motion  of  Sinlaku  for  each  WRF  model  run;  WRE- 
CONV  (a),  WRE-C130  (b),  and  WRE-PSEUDO  (c).  The  red 
line  and  wind  barbs  denote  wind  direction  and  magnitude  in 
increments  of  5  knots. 
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The  times  of  partieular  interest  for  eomparison  of  steering  flows  are  centered 
about  0000  UTC  10  September  because  this  was  when  Sinlaku  began  to  recurve  to  the 
north-northeast  for  a  24  h  period  before  ultimately  turning  back  to  the  northwest  for  its 
approach  to  Taiwan  (Figure  3). 

On  8  September  2008,  the  track  uncertainty  of  Sinlaku  by  the  ECMWF  Ensemble 
Prediction  System  (ECMWF  EPS)  models  was  high  due  in  part  to  the  weak  steering  flow 
of  the  storm,  as  well  as  the  relatively  slow  speed  to  the  north  (Figure  46)  It  was  at  this 
critical  moment  in  the  life  of  Sinlaku  that  would  determine  if  landfall  would  be  in  Taiwan 
or  turn  to  spin  out  harmlessly  over  the  western  North  Pacific. 
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Figure  46.  ECMWF  EPS  at  1200  UTC  8  September.  Pink  lines  represent 
ensemble  member  tracks.  (From  T-PARC  catalog) 

The  predominant  feature  in  all  steering  calculations  was  the  consistent  5-10  kt 
southerly  flow  varying  from  southwest  to  southeast  (Figure  47).  Subtle  details  can  be 
seen  upon  closer  examination  as  it  relates  to  the  different  annulus  sizes.  Flowever, 
regardless  of  the  annulus  size,  it  is  evident  in  WRF-C130  and  WRF -PSEUDO  flow  fields 
prior  to  1200  UTC  9  September  that  a  distinguishable  southwest  component  existed 
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where  none  was  evident  in  WRF-CONV  (Figure  47).  Presumably,  the  faet  that  neither  the 
outflow  jet  nor  the  storm  center  was  well  captured  by  conventional  observations  was 
causes  for  WRF-CONV  to  track  fairly  continually  to  the  west-northwest  without  any 
appreciable  turn  to  the  north. 


lim©  (dd/hh)  fimei  Mrl/hhV 


time  (dd/hht  time  (dd/hh) 


Figure  47.  Average  wind  speed  and  direction  at  300  -  800  hPa  level  with  a 
300-500  km  (left)  and  300-1000  km  (right)  annulus.  The 
abscissa  represents  date  (MM)  and  time  (HH).  The  ordinate  is 
the  magnitude  of  the  steering  flow.  The  green  line  and  wind 
barbs  denote  wind  direction  and  magnitude  in  increments  of  5 

knots. 


When  considering  WRF-PSEUDO  and  WRF-C130,  it’s  important  to  keep  in  mind 
where  the  observations  were  collected.  WRF-PSEUDO  observations  were  collected  well 
north  of  the  storm  center  in  the  upper-level  outflow  region  of  the  storm  while  WRE-C130 
was  collected  at  lower  altitudes  in  close  proximity  of  the  storm  center.  Observational 
studies  have  shown  that  TC  steering  is  often  controlled  by  the  average  air  flow  5-7 
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degrees  from  TC  eenter  (Holland,  1984).  At  this  distanee  the  effect  of  radial  winds  in  the 
steering  flow  calculations  have  little  influence,  thus  isolating  a  more  pure  steering  flow. 

The  flow  in  the  WRF-PSEUDO  300-1000  km  annulus  plots  in  Figure  47  look  to 
be  in  better  agreement  with  the  actual  turn  to  the  north-northeast  that  Sinlaku  took  on  or 
about  10  September  than  does  the  smaller  annulus  of  300-500  km.  The  winds  just  before 
0000  UTC  10  Sept  in  the  300-500  km  annulus  wind  calculations  are  somewhat  opposed. 
In  contrast,  the  winds  in  the  same  time  period  calculated  at  the  larger  annulus  are  more 
parallel  and  complementary  to  the  north-northeast  movement  of  the  storm. 

Ultimately,  the  WRF-C130  experiment  provided  a  more  accurate  TC  track  than 
WRF-CONV  and  WRF-PSFUDO.  However,  it  is  evident  that  a  better  representation  of 
the  overall  environmental  structure;  particularly  where  the  dominant  influence  of  a 
passing  trough  are  concerned,  resulted  in  improved  forecasts.  This  was  evident  in  the 
WRF-C130  and  WRD-PSEUDO  forecast  fields  discussed  earlier. 

By  examining  differences  in  wind  fields  at  200  hPa,  500  hPa,  and  850  hPa  it  is 
evident  that  the  storm  outflow  and  the  midlatitude  jet  were  better  represented  by  the 
WRF-PSEUDO  forecast  while  the  storm  initial  position  was  better  represented  by  WRF- 
C130. 

Comparison  of  the  24-h  WRF-C130  (Figure  48a)  and  WRF-PSEUDO  (Figure 
48b)  850  hPa  forecasts  valid  at  1200  UTC  9  September  indicates  that  easterly  winds  to 
the  northeast  of  the  storm  were  much  stronger  in  WRF-PSEUDO  than  in  the  WRF-C130. 
Furthermore,  this  can  been  seen  in  the  difference  plots  (Figure  48c)  where  a  dipole 
pattern  north  east  of  the  storm  indicates  stronger  easterly  winds  in  the  WRF-PSEUDO 
(dark  blue  shading).  Consequently,  the  WRF-C130  storm  track  proceeded  to  a  more 
northerly  direction  than  WRF-PSEUDO.  The  30-h  forecast  that  verifies  at  1800  UTC  9 
September  forecast  (Figure  49)  indicates  that  the  850  hPa  easterly  winds  in  the  WRF- 
PSEUDO  (Figure  48b)  became  less  in  magnitude  and  shifted  to  more  southeasterly.  This 
was  more  in  agreement  directionally  with  WRF-C130  (Figure  49a),  however  the 
magnitude  of  the  flow  in  WRF-PSEUDO  continued  to  be  stronger  (Figure  49b, c).  In  the 
36-h  forecasts  that  verify  at  0000  UTC  10  September  (Figure  50)  easterly  winds  north- 
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northeast  of  Sinlaku  began  to  weaken  in  the  WRF-PSEUDO  (Figure  50b)  as  the 
subtropieal  ridge  began  to  weaken  and  shift  more  to  the  north  than  in  WRF-C130  (Figure 
50a, e).  This  eontributed  to  the  storm  in  WRF-PSFUDO  to  slow  its  movement  toward  the 
north-northwest. 

The  36-h  forecasts  at  500  hPa  that  verify  at  0000  FITC  10  September  (Figure  51) 
indicate  that  the  midlatitude  cyclone  to  the  north  of  Sinlaku  was  farther  south  in  the 
WRF-C130  (Figure  51a)  forecast  than  in  the  WRF-PSFFIDO  (Figure  51b,c).  This 
southward  extent  of  the  midlatitude  trough  is  consistent  with  the  more  northerly  track  in 
the  WRF-C130  forecast. 

At  200  hPa,  the  24-h  forecast  that  verifies  at  1200  FITC  9  September  indicates 
that  the  outflow  directed  to  the  north  and  east  was  much  stronger  in  the  WRF-PSFFIDO 
(Figure  52b)  than  in  the  WRF-C130  (Figure  52a, c).  This  clearly  demonstrates  the 
potential  effect  of  upper  level  observations  in  the  remote  regions  of  a  TC  as  defined  by 
the  location  of  the  pseudo  dropwindsondes  in  Figure  52b.  By  the  30-h  forecast  that 
verifies  at  1800  FITC  9  September  (Figure  53),  the  merger  of  the  midlatitude  trough  and 
outflow  continued  to  be  stronger  in  the  WRD-PSEFIDO  (Figure  53b)  than  in  the  WRF- 
C130  (Figure  53a,c).  The  stronger  outflow  at  upper  levels  may  be  related  to  the  fact  that 
the  circulation  of  Sinlaku  at  low  levels  is  larger  in  the  WRF-PSEFIDO  forecasts  than  the 
WRF-C130  forecasts.  In  the  36-h  forecasts,  the  pattern  of  stronger  outflow  in  the  WRF- 
PSEFIDO  (Figure  54b)  and  greater  merger  between  outflow  and  the  midlatitude  trough 
than  in  the  WRD-C130  (Figure  54a,c)  continued. 
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Figure  48.  The  24-h  foreeast  winds  (veetors)  and  speed  (m  s-1,  shaded)  at 
850  hPa  and  850  hPa  geopotential  heights  (m,  contours  at  20  m 
intervals)  that  verify  at  1200  UTC  9  September  from  the  (a) 
WRF-C130  and  (b)WRF-PSEUDO,  and  (c)  WRF-C130  minus 
WRF-PSEUDO  forecast  differences. 
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Figure  49.  The  30-h  forecast  winds  (vectors),  and  speed  (m  s-1,  shaded)  at 
850  hPa  and  850  hPa  geopotential  heights  (m,  contours  at  20  m 
intervals)  that  verifies  at  1800  UTC  9  September  from  the  (a) 
WRF-C130,  (b)  WRF-PSEUDO,  and  the  (c)  WRF-C130  minus 
WRF-PSEUDO  forecast  differences  . 
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Figure  50.  The  36-h  forecast  winds  (vectors)  and  speed  (m  s-1,  shaded)  at 
850  hPa  and  850  hPa  geopotential  heights  (m,  contours  at  20  m 
intervals)  that  verifies  at  0000  UTC  10  September  from  the  (a) 
WRF-C130,  (b)  WRF-PSEUDO,  and  (c)  WRF-C130  minus 
WRF-PSEUDO  forecast  differences. 
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Figure  5 1 .  The  36-h  forecast  winds  (vectors)  and  speed  (m  s'\  shaded)  at 
500  hPa  and  500  hPa  geopotential  heights  (m,  contours  at  20  m 
intervals)  that  verify  at  0000  UTC  10  September  from  the  (a) 
WRF-C130,  (b)  WRF-PSEUDO,  and  (c)  WRF-C130  minus 
WRF-PSEUDO  forecast  differences. 


55 


40N 


SON 


20N 


ION 


WRF-C130 


WRF-PSEUDO 


200s-0e-09-1200  -  <4S) 


wknd  sp««d  200  hPa  |m  s''| 

average  geopotanUal  h«lgnt  (m]  200a^>9-09-l200 -(4i) 


w4nd  speed  200  hPa  (m  s') 
average  geopoierttlal  height  |m) 


120E  140e  120E  140E 


200»4>»0»-I200  w4nd  speed  200  hPa  [m  s  'l 

SINLAKU_2  (46)  •  SINLAKU.S  (41)  average  geopoteniial  hoighi  |m] 


120E  140E 


Figure  52.  The  24-h  forecast  winds  (vectors)  and  speed  (m  s'\  shaded)  at 
200  hPa  and  200  hPa  geopotential  heights  (m,  contours  at  20  m 
intervals)  that  verify  at  1200  UTC  9  September  from  the  (a) 
WRF-C130,  (b)  WRF-PSEUDO,  and  (c)  WRF-C130  minus 
WRF-PSEUDO  forecast  differences. 
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Figure  53.  The  30-h  forecast  winds  (vectors)  and  speed  (m  s'\  shaded)  at 
200  hPa  and  200  hPa  geopotential  heights  (m,  contours  at  20  m 
intervals)  that  verify  at  1800  UTC  9  September  from  the  (a) 
WRF-C130,  (b)  WRF-PSEUDO,  and  (c)  WRF-C130  minus 
WRF-PSEUDO  forecast  differences. 
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Figure  54.  The  36-h  foreeast  winds  (vectors)  and  speed  (m  s-1,  shaded)  at 
200  hPa  and  200  hPa  geopotential  heights  (m,  contours  at  20  m 
intervals)  that  verify  at  0000  UTC  10  September  from  the  (a) 
WRF-C130,  (b)  WRF-PSEUDO,  and  (c)  WRF-C130  minus 
WRF-PSEUDO  forecast  differences. 

3,  Storm  Structure  and  Intensity 

The  MSLP  for  each  of  the  experimental  WRF  runs  varied  (Eigure  55).  Unlike  the 
track  forecasts,  WRF-PSEUDO  resulted  in  the  most  accurate  MSLP  forecasts  than  the 
WRE-C130,  or  WRE-CONV.  The  members  of  WRE-CONV  and  WRE-C130  maintained 
a  fairly  constant  spread  until  around  0000  UTC  11  September.  The  WRE-PSEUDO 
members  however  did  not  begin  to  lose  accuracy  until  later  in  the  period  at  1200  UTC  1 1 
September  (Eigure  55c).  It  is  hypothesized  that  the  increased  accuracy  in  storm  intensity 
in  the  WRE-PSEUDO  is  related  to  the  stronger  outflow  from  the  storm.  Perhaps,  the 
addition  of  pseudo  observations  at  upper  levels  in  the  midlatitude  trough  and  outflow 
yielded  better  initial  conditions  of  the  environmental  in  that  key  storm-environment 
interaction. 
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As  with  the  track  forecasts,  by  more  accurately  establishing  the  state  of  the  upper 
level  trough  and  outflow  jet  region  of  Sinlaku  in  WRF-PSEUDO  was  better  able  to 
forecast  the  MSLP  than  WRF-CONV  where  only  conventional  observations  were 
assimilated. 


Figure  55.  Forecasts  of  the  intensity  of  TY  Sinlaku  intensity  for  the  (a) 
WRF-CONV,  (b)  WRF-C130,  and  (c)  WRF-PSEUDO 
forecasts.  The  grey  lines  define  the  ensemble  forecasts  and  the 
orange  line  is  the  ensemble  consensus.  The  black,  red,  green, 
purple,  and  blue  lines  represent  JTWC  best  track,  ECMWF 
YOTC  analysis,  ECMWF  YOTC  forecast,  WRF-ANAF,  and 
WRF-OPS,  respectively.  The  abscissa  is  date  (MM)  and  time 
(HH  UTC).  The  ordinate  is  minimum  sea  level  pressure  (hPa) 
in  5  hPa  increments. 
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V.  CONCLUSION 


The  mission  of  U.  S.  Naval  Meteorology  and  Oeeanography  Command  (METOC) 
has  long  been  eentered  on  providing  the  warfighter  with  the  very  best  environmental 
information  they  need  to  make  the  most  informed  deeision.  To  that  end,  this  thesis  set  out 
to  provide  a  qualitative  example  of  how  a  better  representation  of  the  current  atmosphere 
could  provide  more  accurate  forecasts.  The  western  North  Pacific  relies  heavily  on 
satellite  data  for  atmospheric  observations.  Although  remote  sensing  is  very  accurate  and 
has  been  of  great  service  to  modern  day  meteorologist,  challenges  remain. 

A,  SUMMARY  OF  FINDINGS 

The  goal  of  this  thesis  was  to  build  on  work  done  by  THe  Observing  System 
Research  and  Predictability  Experiment  (THORPEX)  and  more  specifically,  the 
THORPEX-Pacific  Asian  Regional  Campaign  (T-PARC)  that  was  conducted  in  the  fall 
of  2008  (Elsberry  and  Harr  2008).  By  using  the  data  obtained  during  T-PARC,  we  were 
able  to  show  that  with  only  15  additional  upper  level  observations,  strategically  collected 
in  the  outflow  region  of  a  TC  EPS  could  be  benefited. 

Comparisons  were  made  between  three  experimental  WRE  model  forecasts  that 
were  initialized  at  1200  UTC  8  September  and  extended  to  0000  UTC  13  September 
2008.  Eirst,  WRE-CONV  assimilated  only  conventional  observations.  The  second 
forecast  assimilated  conventional  and  WC-130J  dropwindsonde  observations  that  were 
collected  in  Sinlaku  during  T-PARC.  Einally,  WRE-PSEUDO  assimilated  the  same 
conventional  observations  as  WRE-CONV  plus  pseudo  dropwindsonde  atmospheric 
vertical  profiles  derived  from  ECMWF  analysis.  The  pseudo  dropwindsondes  were 
designed  around  a  hypothetic  Global  Hawk  flight  deploying  dropsondes  from  60,000  ft. 
All  WRE  runs  were  initialized  by  ECMWE  operational  forecast  valid  1200  UTC  8 
September,  as  it  would  have  been  the  best  representation  of  the  forecast  if  an  actual  flight 
plan  had  been  performed. 

The  WRE-CONV  forecast  was  less  accurate  than  the  other  two  forecast  models. 
The  WRE  ensemble  consensus  tracked  the  storm  to  the  southern  side  of  Taiwan.  In 
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addition,  the  minimum  sea  level  pressure  was  not  as  low  as  the  other  experimental  runs 
(WRF-C130  and  WRF-PSEUDO).  In  the  WRF-CONV,  the  initial  position  of  the  storm 
was  not  well  characterized.  One  of  the  reasons  for  this  could  be  that  the  conventional 
observations  alone  did  not  capture  the  state  of  the  environment  as  well  as  it  good  have. 

By  examining  WRF-C130  in  light  of  WRF-CONV  a  much  better  initial  placement 
of  the  storm  center  existed  in  the  initial  condition.  This  was  due  in  large  part  by  in  situ 
dropwindsonde  observations  collected  by  the  C-130  reconnaissance  flight.  More 
accurate  initial  placement  served  to  reduce  track  error  over  time  when  compared  to  the 
other  two  forecasts.  Although  the  MSLP  was  closer  to  reality  at  935  hPa,  it  was  still 
higher  than  the  actual  MSLP  of  929  hPa. 

The  initial  placement  continued  to  be  an  issue  with  the  WRF-PSEUDO  forecast. 
Initial  errors  in  the  placement  of  the  storm  center  translated  throughout  the  forecast 
period.  However,  by  providing  a  more  accurate  representation  of  the  upper  atmosphere 
and  the  midlatitude  jet  and  storm  outflow  jet  interaction,  WRF-PSEUDO  was  able  to 
provide  a  more  accurate  MSLP  of  930  hPa  at  0000  UTC  1 1  September.  In  addition,  by 
better  forecasting  the  dominant  steering  trough  to  the  north  of  Sinlaku,  the  ensemble 
consensus  track  resulted  in  a  more  accurate  northern  landfall  in  Taiwan  than  WRF- 
CONV. 

The  results  of  this  study  have  offered  valuable  evidence  into  the  importance  of 
strategically  place  observation  assimilation  into  numerical  weather  prediction.  Although 
the  experimental  pseudo  run,  (WRF-PSEUDO)  in  this  study  fell  short  of  the  track 
accuracy  WRF-C130,  improvements  nevertheless  could  be  seen.  WRF-CONV,  with  its 
sporadic  conventional  observation  DA,  was  improved  upon  by  the  addition  of  a  mere  15 
strategically  procured  observations  in  the  upper  levels  (WRF-PSEUDO).  Also,  WRF- 
PSEUDO  provided  a  much  more  accurate  forecast  of  the  storm  intensity  than  the  other 
run.  This  was  also  owed  to  the  well  captured  state  of  the  upper  levels.  It  is  hypothesized 
that  a  combination  of  observations  over  the  storm  center  and  in  the  environment  would 
provide  the  best  components  of  the  WRF-C130  and  WRF-PSEUDO  forecasts. 
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B,  RECOMMENDATIONS  FOR  FURTHER  STUDY 

Forecasts  of  tropical  cyclone  track  and  intensity  have  long  experienced  large 
uncertainties.  The  sources  of  uncertainty  can  be  numerous.  In  this  study,  a  limited 
experiment  suite  was  conducted  to  assess  the  role  of  observation  distribution  on  forecast 
uncertainty.  In  this  initial  study,  observations  in  the  inner  core  of  the  storm  led  to  the 
most  accurate  track  forecast.  However,  observations  from  the  lower  stratosphere  in  storm 
outflow  and  midlatitude  circulations  lead  to  a  better  intensity  forecast. 

It  is  recommended  that  additional  experiments  be  conducted  to  test  the  value  of 
both  inner  core  and  remote  observations.  All  of  these  observations,  including  those  in  the 
inner  core,  should  be  from  high  altitudes.  Additionally,  additional  experiments  could  be 
conducted  to  test  sensitive  regions  determined  via  methods  used  to  objectively  define 
regions  in  which  there  is  observation  sensitivity  (i.e.,  adjoint  methods,  singular  vectors). 
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